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ABSTRACT 

We develop our model of "punctuated equilibria" for the hot intra-cluster gas 
emitting powerful X-rays. The model considers the gravitational potential wells set 
by the dark matter as they evolve by hierarchical clustering and engulf outer gas; it 
assumes that the gas re-adjusts to a new hydrostatic equilibrium after each merging 
event. Before merging the gas is heated at the virial temperature when bound in 
subclusters; at early z it is preheated by supernova activity following star formation. 

In detail, we compute analytically the following steps: the dynamic histories of 
dark matter halos with their merging events; the associated infall of gas into a halo, 
with compressions and shocks estabilishing the conditions at the cluster boundary; 
the updated disposition of the gas in the potential well matching such conditions; the 
statistical convolution of observable quantities over the merging histories. 

For the individual objects from groups to clusters, the model yields profiles of 
density and surface brightness with no free parameters; in particular, the so-called /3 
parameter is itself an outcome of the model, and the polytropic index 7 is internally 
constrained to a narrow range. We obtain declining temperature profiles, and profiles 
for the density and for the surface brightness shallower in groups compared with 
clusters; our model groups also contain a lower baryonic fraction on average, but with 
a scatter considerably larger. 

We present various key quantities over the whole range from groups to clusters. 
In particular, we predict in different cosmologies the statistical correlation L — T of 
luminosity with temperature; similarly, we derive the correlation Rx — T for the size of 
the X-ray emitting region. The intrinsic scatter in both correlations is also predicted. 
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1. INTRODUCTION 

Groups and clusters of galaxies constitute large, nearly virialized condensations. Within their 
virial radii R ranging from 1/2 to about 2 l^pc the density contrasts attain or exceed Sp^pu ^ 
relative to the background, and the corresponding masses M range from some 10^^ to 10^^ Mq, 
mainly in dark matter (DM). 

These structures contain a large baryonic fraction / 2± 0.1 in the form of a hot intra-cluster 
plasma (ICP) at temperatures kT ^ GM mn/^OR ~ 0.5 — 10 keV {niH is the proton mass) with 
particle densities up to n ~ 10~^ cm~^. The ICP emits by thin thermal bremsstrahlung X-ray 
luminosities L oc n ranging from 10^^ to 10^^ erg/s. 

Recent observations of individual objects resolved in energy and angle with the ASCA and 
with the SAX satellites indicate radial temperature profiles declining outwards (Markcvitch, 
Sarazin &: Henriksen 1997; Molendi, private communication; see also Hughes, Gorenstein &; 
Fabricant 1988). In several cases the thermal structure is complicated by hot spots (Honda et 
al. 1997; Markevitch et al. 1998). As for the statistical aspects, a steep correlation close to the 
overall form L oc is known to hold for local clusters, but with a substantial scatter (Edge &; 
Stewart 1991; Mushotzky 1994). Recently the observations have sampled higher redshifts out 
to 2: ~ 0.5 (Tsuru et al. 1996; Mushotzky &: Scharf 1997) finding little significant evolution. 
Low-temperature, local systems have been also sampled, finding there indications of a slope 
steeper yet (Ponman et al. 1996). For kT > 5 keV a flattening toward L ocT^ has been detected 
by Allen & Fabian (1997). 

In the near future the AXAF mission will substantially improve the space-resolved 
spectroscopy of many individual clusters, and shortly after the mission XMM will greatly enlarge 
the statistics of the L — T correlation. Corresponding upgrades are called for in the theoretical 
understanding of the complex astrophysics concerning both the DM and the ICP over the full 
range from groups to clusters. 

The force approach uses numerical computations for both the DM and the ICP, striving 
for wide dynamic range and complete hydrodynamics. The N-body simulations have shown, in 
accord with the hierarchical clustering picture (sec Peebles 1993), the evolution of the DM halos 
to occur largely through a sequence of merging and accretion events which involve generally 
smaller partners down to nearly diffuse matter (see, e.g., Tormen, Bouchet & White 1997), and 
are correlated with the surrounding large scale structures (Colberg et al. 1998) 

As for the ICP, pioneering work by Schindler &; Miiller (1993) taken up by Roettiger, Stone 
&; Mushotzky (1998) used 3D Eulerian codes with adaptive mesh and advanced shock-capturing 
methods to study how the large merging events of the DM halos affect the ICP component. 
The outcomes show how such events produce anisotropic shocks and non-uniform compressions, 
resulting in a complex thermal structure lasting a few Gyrs. 

At the other extreme, a sequence of radial, highly resolved Eulerian computations (progressing 
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from Perrenod 1980 to Takizawa &; Mineshige 1998a) have shown that isotropic accretion of 
smooth gas also causes a strong and slowly expanding shock, which remains close to the (growing) 
virial radius for some dynamical times. 

Most recently, state-of-the-art N-body codes coupled with advanced hydro (Bryan Sz Norman 
1998, Ghcllcr ct al. 1998) have been run on supercomputers, aiming at resolutions below 100 kpc in 
rich clusters as necessary for deriving reliable luminosities. But such simulations are hard-pressed 
in implementing at the same time the full physics of the ICP. In fact, preheating at temperatures 
~ 0.5 keV is expected from stellar formation and evolution to supernovae (see Renzini 1997, and 
references therein); this is particularly relevant for the shallower potentials of groups where T is 
close to 0.5 keV. Inclusion of the preheating in the numerical work is technically taxing, as it 
involves cooling, star formation and energy feedbacks resolved down to subgalactic scales; but 
in its absence the simulations produce a correlation of the form L oc at all temperatures, at 
variance with the data. Suginohara Sz Ostriker (1998) stress how delicate may become at high 
resolutions the balance of cooling and feedbacks, and how difficult becomes reproducing ICP cores 
as observed. Thus for now and for some time to come it will be hard to combine into a realistic 
numerical picture wide dynamic range from galaxies to large scale structures, stellar preheating 
and large statistics. 

The state of the numerical approach and the challenge from the data motivates us to present 
here an analytic model which includes, though in a simplified form, the physical processes outlined 
above. We describe the cluster history as a sequence of punctuated equilibria (PE). That is to 
say, we envisage such history as a sequence of hierarchical merging episodes of the DM halos 
which we compute analytically (with its variance) in the framework of the standard hierarchical 
clustering, specifically using the so-called "extended Press &; Schechter theory" (Bond et al. 1991; 
Bower 1991; Lacey & Cole 1993). We stress that such episodes cause in the gas shocks of various 
strengths depending on the mass ratio of the merging subclusters, ranging from nearly adiabatic 
compressions for comparable clumps up to shocks with high Mach numbers in the accretion of 
loose gas. Our point is that the most effective such shocks and compressions overlap to provide the 
boundary conditions for the new hydrostatic equilibrium to which the ICP is assumed to re-adjust. 

The PE model as presented here takes up our previous work (Cavaliere, Menci Sz Tozzi 
1997, 1998; hereafter CMT97, CMT98), but differs in that the histories of the DM halos are 
now computed analytically rather than based on Monte Carlo simulations. This goes beyond the 
technical aspect since it allows us to explore efficiently the dependences of the density, temperature 
and luminosity on the parameters of the clusters and on the cosmology. In the same vein, the 
present approach allows us to quantify the connection between slope and scatter of the L — T 
correlation and the cosmological scenario. In addition, the parameters of the ICP thermal state 
are now fixed or bounded in terms of constraints internal to the model, and new predictions are 
presented. 

In Sect. 2 we describe the PE model and our computational steps: in Sect. 2.1 we recall the 
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statistical formalism for hierarchically merging DM halos; in Sect. 2.2, we derive from the mass 
ratios involved in each merging episode the boundary conditions for the plasma equilibrium; in 
Sect. 2.3 we compute from such boundary conditions the TCP equilibrium; in Sect. 2.4 we derive 
the statistics of L and of the size Rx using the formalism of Sect. 2.1. In Sect. 3 we give, and 
compare with the observations, the model results for the profiles n(r) and T(r), for the surface 
brightness S(r), for the relations M — T and f — T, and for the correlations Rx — T and L — T. 
The final Sect. 4 is devoted to discussion and conclusions. 



2. THE PUNCTUATED EQUILIBRIA MODEL 

The X-ray bolometric luminosity of a cluster is given in its basic dependences by 

Loc n n^{r)T^/'^{r)Sr . (1) 



Here T(r) is temperature in the plasma and r2 is the cluster boundary, that we take to be close 
to the virial radius R oc M^/^p"^/^, where p{z) oc (1 + z)^ is the DM density in the cluster, 
proportional to the average cosmic DM density Pu{z) at formation. 

It will be convenient to separate the internal profiles n(r) and T{r) from their boundary 
conditions at r2. As to the latter, the infalling gas is expected to become supersonic near r2 (see, 
e.g., Pcrrcnod 1980; Takizawa & Mincshigc 1997) so that a shock front will form there. The 
conservations across the shock of mass, energy and stresses yield the Rankine-Hugoniot conditions, 
i.e., the temperature and density jumps from the outer values Ti and n\ to T2 and n2 just interior 
to r2 (spelled out in Sect. 2.2). Then the luminosity may be rewritten in the form 

Lo^rlniT^ f^'.[^]'m]''\ (2) 

JQ L n2 -I L i2 J 

where x = rjr^- 

Note that the values n2 and T2 at the boundary are not uniquely determined by the cluster 
mass M; rather, they are related to the outer values n\ and T\ by the named shock conditions. 
In turn, n\ is fixed by n\ oc fuPu/f^H, in terms of the universal baryonic fraction /„; whereas 
Ti is determined only statistically, through the diverse merging histories ending up in the mass 
M. Specifically, as explained in detail in Sects. 2.1 and 2.2, in each merging episode Ti takes 
on the values appropriate to the other merging partner, constituted by subclumps or even by 
smooth gas. In sum, a given dark mass M admits a set of ICP equilibrium states characterized 
by different boundary conditions, each corresponding to a different realization of the dynamical 
merging history. It is the convolution over such set which provides the average values of L and 
Rx, and their scatter. 



So the development of our PE model proceeds along the following steps: 
i) we first give the statistics of the current DM halo of mass M and of the merging clumps AM; 
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ii) we compute the shock strength relating at the boundary the inner values T2 and n2 to the 

exterior ones Ti and ni, as a function of M and AM; 

iii) from such boundary values, we compute the internal profiles T{x)/T2 and n{x)/n2 for the 
post-merging hydrostatic equilibrium, involving the cluster potential and hence M; 

iv) we convolve the results of steps ii) and iii) with the statistics i). 

Below we describe these steps in turn. 



2.1. Histories of the Dark Matter Halos 



Here we recall the basic merging probabilities provided by the "extended Press & Schechter 
theory", see Lacey & Cole (1993; 1994). This is based on the dark halos formed by hierarchical 
merging of smaller structures. 

The halo mass distribution at the cosmic time t is given by the standard Press & Schechter 
(1974) formula 



N{m,t) 



2 5c{t)pu 
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(3) 



where the masses m = M/Mg are normalized to the current value Mq = 0.6 10^^ Q,o h"^ Mq (i.e., 
to the mass enclosed within a sphere of 8h^^ Mpc), and 6c{t) = 6coD{t) is a critical threshold for 
the collapse and virialization of the primordial density perturbations. The local value 6co depends 
weakly on the cosmological parameters, while the growth factor D{t) sensitively depends on them. 
The mass variance a{m) is computed in terms of the perturbation spectrum; for definitness, we use 
the CDM spectra given and discussed by White, Viana &; Liddle (1996). For f2 = 1 we adopt the 
primordial "tilted" index Up = 0.8; for 0,o < 1 we adopt rip = 1. The associated normalizations are 
taken from the COBE/DMR data (Gorski ct al. 1998), and expressed in terms of the amplitude 
(Tg at the relevant scale of 8h^^ Mpc (sec Bunn &: White 1997). 

Corresponding to eq. (3), the probability distribution that a given mass m at time to has a 
progenitor of mass m' at time ti < to reads 
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(4) 



where cr' is the mass variance at the scale m'. On the other hand, at a given time t a progenitor 
m' increases its mass by a merging event with a clump of mass Am (producing a cluster with 
mass m = m' + Am), with the probability distribution per unit time given by 

d'^p{m' ^ m' + Am) _ f2\i/2 dln{Sc) dlna 



dAmdt 



(- 
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(5) 
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We have compared the analytical probabilities above with the results of the Monte Carlo code 

developed by P. Tozzi to simulate the hierarchical merging history of halos, based on the excursion 
set approach of Bond et al. (1991). A realization from the Monte Carlo simulations is shown 
in fig. la as an illustration of the basic process of DM halo growth. To show the agreement 
of the two approaches for a relevant quantity, we plot in fig. lb the probability distribution of 
progenitors of mass m' at different z that end up in a given mass m at z = 0, computed from the 
Monte Carlo and from the eqs. above. In fig. Ic we show as a function of m, the fraction of objects 
which, during the last 2 Gyrs, accreted mass in events involving comparable clumps (specifically, 
those with mass ratio 1/2.5), and in events involving very unequal clumps (with mass ratio 1/10); 
during that interval, more than 60 % of the clusters with M > 10^^ Mq will have merged with 
clumps smaller than M/10. 



2.2. Boundary Conditions 

The pre-shock temperature in a merging event is that of the infalling gas. If the latter is 
contained in a sufficiently deep potential well, Ti is the virial temperature Ti^ oc i\.mlr of the 
secondary merging partner; on using r oc {^t^mj pf-l^ this writes 

kTxv = 4.5 (Am)2/3 (p/po)V3 keV, (6) 

where the numerical coefficient is taken from Hjorth, Oukbir &: van Kampen (1998). Where 
necessary, the z-dependence of p/po = (1 + z)^ is converted to ^-dependence following the standard 
FRW cosmologies. 

But an independent lower bound /cTi* « 0.5 keV is provided by preheating of diffuse external 
gas, due to feedback energy inputs following star formation and evolution all the way to supernovae 
(David et al. 1995; Renzini 1997). We recall that preheating temperatures in excess of 0.1 keV 
are believed to constitute essential complements to the hierarchical clustering picture to prevent 
the "cooling catastrophe" from occurring, see White & Rees (1978); Blanchard, Vails Gabaud & 
Mamon (1992). In point of fact, Henriksen & White (1996) find from X-rays evidence for diffuse 
gas at 0.5 — 1 keV in the outer regions of a number of clusters. So in the following the actual value 
of Ti will be 

n = max[ny,Tu]. (7) 

Given Ti, the boundary conditions for the ICP in the cluster is set by the strength of the 
shocks separating the inner from the infalling gas. We report here from CMT98 the explicit 
expression of the post-shock temperature T2 for three degrees of freedom and for a nearly 
hydrostatic post-shock condition with V2 « vi, assuming the shock velocity to match the growth 
rate of the virial radius R(t): 



4 10 2O(l + 0:Ti)2 



(8) 



- 7- 



Here e = 15kTi/4fimHVi and n is the average molecular weight; the inflow velocity vi is set 
by the potential drop across the region of nearly free fall, to read vi ~ \J —^21^^11 in terms of 
the potential (^2 at r2- For a "cold inflow" with e << 1 the shock is strong, and the expression 
simplifies to kT2 ~ + 2>kTi/2. Instead, for e 1 the shock is weak, and T2 ^ Ti is 

recovered as expected. Note that T2 depends through both Ti and v\ on the mass Am of the 
clump being accreted. 

^From T2 and Ti, the density jump at the boundary n2/ni is found to read (see CMT97) 

ni V Ts/ L V T2^ ^ TsJ ^ ' 

It is seen that the density jump takes on the limiting value n2/ni = 4 for very strong shocks, while 
the adiabatic approximation n2/ni « 1 + 2>{T2 — Ti)/2Ti is recovered for weak shocks. 



2.3. Hydrostatic Equilibrium 



We adopt the polytropic temperature description T{x)/T2 = [n{x)/n2\^ ^, with the index 7 
in the range l<7<5/3to begin with. In terms of the virial temperature (see Sarazin 1988), 
eq. (2) writes L oc {ni/pf (n2/ni)2 T^^ p^^Ta/T^)^/^ [n(r)/n2]^'*"^'^~^^/^ where the bar denotes 
the integration over the emitting volume < rf, and p is the average DM density in the cluster, 
proportional to Pu and so to ni (see under eq. 2). The radius r2 may be rewritten in terms of the 
temperature T^, oc m/r2 oc pr^. We finally obtain 



1/2. 



[n{r)/n2] 



2+(7-l)/2 



(10) 



The underlying assumption is that after a merging event the cluster re-adjusts to a hydrostatic 
equilibrium with boundary conditions n2, T2 corresponding to its dynamical history (see eq. 
6-9). Actually, this requires sound propagation times shorter than the dynamical timescale taken 
anyway by the DM to relax to its own steady configuration; the condition is seen to be satisfied, 
though marginally, for all merging events except for the rare ones involving comparable clumps. 
As we discuss in detail in the concluding §4, the observations and the hydrodynamical N-body 
simulations concur in supporting not only the hydrostatic equilibrium approximation for the 
relevant merging events, but also its parametrization with a polytropic equation of state. 

The ratio n{x)/n2 is obtained from the hydrostatic equilibrium dP/rriH ndr = —GM{< 
r)/r'^ = —d<p/dr with the polytropic pressure P(r) = /cT2 n2 [n{r)/n2]^- This yields (see Cavaliere 
&; Fusco Femiano 1978; Sarazin 1988, and bibliography therein) the profiles 
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T(r) 



T2 



L 7 
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(11) 
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where = (l)/iJ,mH (^2 is the potential normalized to the associated one-dimensional DM velocity 
dispersion at r2- The ICP disposition in eq. (11) relative to the DM depends on the parameter 

/3 = iimH(T2/kT2 , (12) 

and is further modulated by the second parameter 7, to yield as the latter increases flatter profiles 
n(r) and steeper T(r). In our PE /3 is by T2 given by eq. (8), and considering the statistics of T2 
we obtain the results discussed in §3.1; the other parameter 7 will be bounded as also discussed 
there. 

We shall focus on the "universal" forms of 0(r) and a{r) given by Navarro, Frenk &; White 
(1997). When relevant, we will discuss also results for the (simplified) King potential (see Sarazin 
1988; see also Adami et al. 1998) where the DM itself has a core, and somewhat fatter ICP cores 
obtain. We shall also discuss the steeper cusp found by Moore et al. (1997) in highly resolved, 
CDM simulations; correspondigly, we still obtain a core-like ICP distribution, albeit slightly 
slimmer. Actually, in hydrostatic equilibrium a DM cusp flatter than p{r) oc - corresponding 
to a gravitational force flatter then - implies at the centre a finite ICP density ric but a high 
derivative, which however at observable resolutions is flattened by a modest increase of 7. 



2.4. Statistics 

Our purpose is to compute the average value of L and its dispersion, associated with a given 
cluster mass m. We re-iterate from Sect. 2 that the diverse merging histories ending up in such a 

mass give rise to a set of equilibrium states characterized by different values of T2, n2- These are 
related by eqs. (8) and (9) to the values of Ti associated with the the clump of mass Am incoming 
onto a cluster progenitor. So to meet our purpose we must sum over the shocks produced at a 
time t' < t in all possible progenitors m' (weighting with their number) by the accreted clumps 
Am (weighting with their merging rate); flnally, we integrate over times t' from an effective lower 
limit t - At. 

The average L is then given by 

(L) = Q / dt' / dm' / dAm-^(m',l/\m,t)—^ — ^L; (13) 

^ ' Jt-At Jo Jo dm'^ ' I w dAmdt' ^ ' 

and the variance is given by 

,A.2^ ^ /■* ,//■"",/ /■'""'"' ,A df , . ,, , d^pim' ^ m' + Am) / ^ ,,,\2 

(14) 

Higher moments - if needed in case of non-Poissonian statistics - are given by similar expressions; 
the full distribution of L requires aimed computations or simulations, as noted by CMT98. In 
the integrals, the luminosity L [r2(m'), ri(Am')] depends on m' and Am' through the boundary 
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conditions discussed in Sect. 2.2. The compounded probability distribution in eqs. (13) and 

(14) has been normahzed to 1 (we do not write down the normahzation factor Q for the sake of 
simphcity) . The effective lower limit for the integration over masses is set as follows. 

The merging events relevant to (L) and to (AL^) after eqs. (13) and (14) are those lasting 
enough as to overlap with similar ones. Since At oc Arn/^i pr^, the above condition results in 
an effective lower limit for the masses entering eq. (13) and (14); physically, lumps with masses 
smaller than such limit yield a small mass flux and so produce shocks which dissipate before new 
clumps income. We pinned down the minimum Am numerically, by looking at the saturation of 
(L) for increasing values of At. This occurs at At k. 0.7 td which corresponds to a lower limit 
about m/20 for Am on using oc Am^/^. 

Note that the above procedure acts like an effective mass weight. Heuristically, this may 
be seen with the Am and t integrations interchanged; then the lower limit contains At (Am) 
which must be convolved with the distribution of Am. But for At > 0.7 1^ the resulting average 
saturates; so we have written the t-integration as the outmost integral, having adopted the lower 
limit for Am said above. 

Thus very small accreted lumps do not affect our average {L), due to the small associated 
mass accretion rate. On the other hand, merging events involving comparable partners (though 
contributing ~ 1/2 of the total mass) affect the overall (L) only marginally; in fact, such events 
not only are few (< 10%), but also they involve lumps with temperatures Ti comparable to T2, 
and so produce a compression factor 712/^1 ~ 1 (eq. 9). The major contribution to (L) is by far 
(~ 90%) provided by intermediate merging lumps, which yield an integrated contribution ~ 1/2 to 
the mass, but dominate the number of events and produce large compressions (^2/71-1)^ oc L. For 
such events, the isotropic hydrostatic equilibrium for the ICP is physically motivated and robust 
(see §2.3 and §4), substantiating our step-by-step rendition of the hydrodynamics. 

3. RESULTS 

Here we present various results from the PE model, and compare them with observations. 
To this aim, we shall express our results in terms of the observed emission weighted temperature, 
which we denote simply by T. Moreover, the contribution of relevant emission lines to L (from 
updates of Raymond & Smith 1977) has been added to the bremmstrahlung emission underlying 
the simple scaling in eq. 10. 

3.1. Profiles 

Our reference cluster will have a mass m, and DM potential 0(r) as said in Sect. 2.3. The 
density and temperature profiles are given by eq. (11), and are to match the shock boundary 
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conditions at the position — R- The key quantity is the parameter (3 defined by eq. (12); its 

average value and scatter are predicted by the PE model (iising convolutions analogous to eqs. 
13-14) to be as shown in fig. 2. Note that /3(T) grows slowly with the temperature; we obtain 
values ranging from about (5 = 0.5 at the group scale to /3 « 0.9 for rich clusters, see fig. 2. These 
values imply that in our model the ICP profiles are smoother and more extended than the DM's, 
an effect becoming more pronounced in downscaling from clusters to groups. 

In fig. 3 we show the baryonic fraction /2 (integrated out to the shock) as a function of T and 
for different values of 7. The polytropic index 7 > 1 describes the equation of state for the ICP. An 
upper bound to it arises if the overall thermal energy of the ICP is not to exceed its gravitational 
energy, with only minor contributions from other energy sources, like radiosource heating or energy 
transfer from DM to ICP, as discussed in Sect. 4. The thermal and the gravitational energy are 
computed using the profiles in eq. (11), and their ratio is given in fig. 4 to show that the upper 
bound 7 < 1.3 holds. 

In fig. 5 we show temperature profiles T(r) for different values of 7, in terms of the normalized 
coordinate x = r/R. It turns out that observations by Markevitch et al. (1997) are consistent with 
the r(r) predicted when 7 = 1.2 it 0.1, in our allowed range. Hereafter we shall focus on 7 = 1.2. 

In fig. 6 we show the density profiles n(r) for two local clusters with different temperatures; 
the associated surface brightness E(r) is shown in fig. 7 along with representative data. It can be 
seen that groups have flatter E(r) than rich clusters, an outcome persisting when the King or the 
potentials by Moore at al. 1997 are used. 

3.2. Correlations 

We show first in fig. 8 the M — T relation in view of its important role. Note that, given the 
mass function, our flattening at low temperatures translates into a steepening of the corresponding 
temperature function; such an effect has been also noted by Balogh, Babul & Patton (1998). 

The L — T correlation is given by the double convolution (13), and likewise for AL after eq. 
(14). The results are shown in flg. 9 for a tilted CDM spectrum of perturbations in the critical 
universe. For the reasons discussed in CMT98, the normalization has been best-fitted on the data. 

As stressed in our previous works (CMT97, CMT98), the correlation we predict and show in 
fig. 9 is not a single power-law; it starts as L oc for very rich clusters with high T, but bends 
down with decreasing T, due to the threshold effect of the preheating temperature kTi ^ 0.5 keV. 
Note that our correction to 0.3 solar metallicity tends to increase, if anything, the luminosities 
and to flatten the slope at low T. 

A convenient fltting formula for the predicted L — T correlation (precise to better than 10% 
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for T > 2ri) is as follows: 

L = aLT'^+'^Hp/pof' (15) 
aL oc 3.8J^O-=^(l + zf2Vf^° + (l-Oo)e-°-^(^+^) 
aL = i.i2(l + 2)-0-2e-0-25(^-^i)/f^r(i+^)°-\ 

where the luminosity is expressed in units of 10^^ erg/s and the temperature in keV. The 
z-dependence of results from the interplay of the following effects: i) the evolution of the 
Navarro et al. (1997) potential; ii) the abundance of clusters with given T; iii) the evolution of 
the progenitor probability distributions in eqs. (4) and (5). Such effects are small, and moreover 
they balance out very nearly, leaving the basic z-dependence [p{z) / po]^/"^ . In turn, the latter 
dependence goes as (1 + z)^'^~^, considering (CMT98) the evolution of the overdensities in the 
large scale structures - filamentary or sheet-like - hosts to the clusters. 

At temperatures substantially larger than the threshold kTi ~ 0.5 keV the intrinsic, dynamic 
dispersion grows with T, but the relative AL/L stays nearly constant around 25 % (20"), as shown 
by fig. 9. We have deliberately chosen to keep this figure simple and not to include the conceivable 
spread of Ti already represented by CMT98 in their fig. 2; the effect of such a spread is to widen 
the dispersion below 1 keV adding another, large component to the intrinsic dynamic variance. 

As our analytical approach allows us to span a wide range of cosmologies /cosmogonies, we 
show in fig. 10a the dependence of (L) on at two temperatures and at the current epoch. 
It is seen that {L) increases with Qq increasing; this is because the underlying strength of the 
current shocks grows on average as the merging rate (moderately) increases on approaching the 
critical cosmology, see Lacey & Cole (1993). A similar behaviour is followed by the corresponding 
dispersion AL, see fig. 10b. 

Similarly, we derive a correlation with T of the effective size of the X-ray emission. If 
(following Evrard & Mohr 1997) we define Rx in terms of the isophote corresponding 1.9 10""^ 
ct/s arcmin^ in the ROSAT band (consistent with our normalization for L), we find the {Rx) — T 
correlation shown in fig. 11. We find Rx oc T in the range of clusters, with a steepening at the 
group scale and a flattening at large temperatures. A corresponding fitting formula is as follows: 

Rx = aflT0-5+"«(p/p„)-V2 (16) 

UR = 0.6n°-i(l + z)-i-5+°-=^(i-^o)/(i+^) 

aR = 0.5V'(l + ^)°-'e-°-=^'^^"^^^''°"/^'+'^" > 

where Rx is expressed in Mpc and T in keV. 

4. CONCLUSIONS AND DISCUSSION 

This paper is based on hierarchical clustering; group and cluster formation is envisaged in 
terms of DM potential wells evolving hierarchically, and engulfing outer baryons by accretion of 
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smooth gas or by merging with other clumps. We consider the diffuse baryonic component to 

increase as the deepening wells overcome the external gas energy provided by preheating stars 
or by virialized subclumps. After a merging episode, the TCP in the wells falls back to a new, 
approximate hydrostatic equilibrium (hence the name "punctuated equilibria"). 

We have modeled this complex astrophysics using an analytic approach based on the standard 
hierarchical clustering and comprising three main steps: the hydrostatic equilibrium for the ICP 
is computed for a given boundary condition; the latter is derived from the effects on the ICP of 
the dynamical evolution of the DM halos; the intrinsic stochastic character of such evolution is 
accounted for on convolving with the statistics of the DM merging histories. All such aspects 
are treated in a fashion which is necessarily simplified; however, the resulting model - whose 
parameters are all internally or physically constrained - proves to be gratifyingly efficient in 
explaining and predicting a variety of observations. We shall discuss such aspects in turn. 

As for the ICP equilibrium, we note that the condition R/cg < (sound crossing time 
shorter than the dynamical time) is weakly satisfied. However, this is enough to ensure ICP 
equilibrium during the intervals (at least 2/3 of the total time, see Tormen et al. 1997) when the 
DM halos themselves are in approximate dynamical equilibrium. In point of fact, ever since Jones 
&; Forman (1984) to the recent Cirimele, Nesci Sz Trevese (1997), it has been recognized that 
hydrostatic equilibrium provides for many clusters a fitting description of the averaged profiles of 
surface brightness in X-rays (except for the central region when cooling fiows set in). Even the 
high-resolution observations provided by ASCA (Markevitch 1998) and those being derived from 
SAX data can be accounted for in terms of average profiles. Conspicuous hot spots do occur, but 
only in a minority of sources, and then in correlation with other signs of ongoing major dynamical 
events, as discussed next. 

This body of evidence supports the case that the gas stays close to hydrostatic equilibrium, or 
falls back to it after a short transient from the merging event, except for the rare major episodes 
involving a partner of comparable mass. The limits to the above picture may be defined with the 
help of aimed hydro dynamical simulations. A quantitative account of how much and how long 
cluster collisions displace the ICP out of equilibrium can be found, e.g., in the N-body experiment 
of Roettiger et al. (1998) for the case of a merging event with a mass ratio of 1/2.5. Even for such 
ratio (already a rare event in the hierarchical clustering picture) the simulation shows that some 2 
Gyr after the event hot spots and space variations of the luminosity are reduced to under 20%. 

So a sequence of hydrostatic equilibria of the ICP is physically motivated for all merging 
events except for those involving comparable clumps (a mass ratio larger than ~ 1/4). However 

these sum up to less than 10% in the number; in addition, these events yield a shock compression 
factor n2/ni 1 (see eq. 9), with an overall contribution to (L) less than 10 %. This is actually 
the precision level of our model. 

Note that these considerations also support the use of the polytropic relation T ~ n'^~^. In 
fact, when equilibrium holds, a macroscopic 7 = dlnp/dlnn may be always defined in principle to 
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describe the ICP state; the question concerning whether this is constant on scales ;^ 0.1 Mpc can 

be probed with observations. In fact, as we discuss later on in this Section, observed temperature 
and surface brightness profiles agree well with those predicted by the polytropic equation of state 
with 7 ~ 1.2. 

The hydrostatic equilibrium is described by a first order differential equation (sec §2.3) 
requiring one condition at the boundary. Physically, this is provided by the place where a shock 
converts most of the kinetic and gravitational energy of the inflowing colder gas into thermal 
energy, as it must occur for the accreted ICP to be contained in the well. The boundary condition 
may be referred to in terms of the stress balance P2 = Pi + nimnVi, one conservation law 
contributing to eqs. (8) and (9). If it were somehow possible to shut off the r.h.s. completely, 
the intracluster medium would expand (the more so the closer is its state to isothermal, and the 
shallower is the potential at the shock position; the density would decrease everywhere including 
the centre and L oc would quench considerably over a sound crossing time. However, we shall 
argue next that the dynamic stress acts steadily. 

In closer detail, the shock jump conditions are set in terms of T2/T1, basically the height of 
the current potential well (see eq. 8) compared with the thermal energy of the infalling gas. The 
latter is initially due to stellar preheating (of nuclear origin); then it is increased to the virial value 
(of gravitational origin) when the accreted gas is bound in DM subclumps. So the preheating sets 
an effective threshold kTi ~ 0.5 keV to gas inclusion, which breaks the self-similar correlation 
L (X not only in its vicinity but also up to a few keV. In our model, this occurs through the 
specific dependence of n2/rai on T2/T1 at the cluster boundary holding for spherical shocks, strong 
or weak. This is a fair representation for the conditions prevailing when the cluster growth occurs 
by nearly isotropic accretion of smooth gas or of many small clumps, as shown by a sequence 
of spherical hydrodinamical simulations up to the recent one by Takizawa & Mineshige (1997). 
Surely, this representation looks as a rather crude approximation to the aftermaths (lasting up 
to 2 Gyr) of major merging events, as those simulated in detail by Roettiger et al. (1998). But 
in point of fact, our average quantities, their scatter and the profiles agree with the data over the 
whole range of T as we stress next, while an explanation in point of principle is offered thereafter. 

For example, the equilibrium parameter /3{T) is set by the boundary conditions in terms 
of shock strengths to values which increase from about 0.5 at the group scale to about 0.9 for 
rich clusters, see fig. 2. Correspondingly, the surface brightness profiles in groups - beyond the 
generally larger observational noise - ought to be flatter than in rich clusters (see figs. 6, 7). In 
fact, similar values are obtained from fits to the brightness S(r) observed in groups and in rich 
clusters ever since Kriss, Cioffi &; Canizares 1983 and Jones &: Forman 1984. More recently, a 
similar trend has been found from spectroscopic measurements of (3 by Edge &; Stewart (1991), by 
Henriksen et al. (1996), and by Girardi et al. (1998). We add that - as another straightforward 
consequence of the threshold Ti - our model groups differ from rich clusters also for their lower, 
average baryonic fraction (see fig. 3) in accord with the average values inferred by Reichart et al. 
(1998) from observations, however noisy. 
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The other parameter of the ICP equihbrium, namely the polytropic index 7, is constrained 
to the rather narrow range 1 < 7 < 1.3. The upper bound holds if the ICP thermal energy 
content is not to exceed its gravitational energy in the DM potential (see figs. 4); it may be 
extended to 1.4 only if kinetic energy transfer from residual clumps in DM to ICP contributes 
more than 20%, the limit set by Kravtsov & Klypin (1998). Values of 7 around 1.2 imply the 
entropy distribution S{r) oc log[T{r)/n'^~^{r)] to have a neat central minimum, in accord with 
the notion of dominant entropy deposition by shocks in the outer regions (see David, Jones & 
Forman 1996; Bower 1997) against the central contribution deposited by supernovae. In fact, 
specific calculations based on entropy production at the shock show that 7 ~ 1.2 holds with little 
variations from clusters to groups (Tozzi & Norman 1999, in preparation). Values of 7 ~ 1.2 turn 
out to yield temperature gradients (see fig. 5) consistent not only with the results from advanced 
simulations of rich clusters (Bryan & Norman 1998), but also with aimed recent observations by 
Markevitch et al. (1997) (see also Fusco Femiano & Hughes 1994). Preliminary data from SAX 
(S. Molendi, private communication) indicate in some cases a somewhat flatter central gradient, 
but still consistent with the range of 7 given above. That such gradients cannot be realistically 
traced back to imperfect thermalization of the electrons has been argued by Ettori & Fabian 
(1998) (for a discussion see also Takizawa & Mineshige 1998b). Note that the arguments may 
be reversed, opening an interesting perspective to gauge the baryon thermal history and its link 
with the DM dynamics, see fig. 6; in fact, very steep DM cusps would require larger values of 7 
to fit the core-like shape of S(r), but this in turn would produce steep profiles of T{r) and imply 
additional central inputs of energy and entropy, leading to a strongly bent L — T relation. 

As for the statistical aspects governing average correlations and their scatter, these are 
derived from convolution of the boundary conditions over the merging histories. We have already 
predicted and discussed the L — T correlation in our previous papers (CMT97; CMT98); here 
we only note, and illustrate in fig. 9, the points discussed by Markevitch (1998), Allen &; Fabian 
(1998), namely, that once the effect of large cooling flows is removed or accounted for, the average 
correlation is flattened to a slope around 2.5 and the scatter is reduced down to 13 % (Icr), both 
in good agreement with the intrinsic, dynamic scatter from the model. 

Here instead we expand on the Rx — T correlation, see fig. 11. That our model predicts 
Rx oc T in the range of rich clusters with a steepening at the group temperatures, is due both to 
the non selfsimilar form of the shock strength as a function of T, and to the shape of P{T) discussed 
above. The results agree with the data, whilst all the self-similar computations (including the 
simulations without preheating discussed by Evrard &; Mohr 1997) yield Rx oc rO-5-o.7 rpj^g 
dispersion we find from the average over the merging histories also compares well with the existing 
data. 

To understand such overall effectiveness of the model one has to consider two features of the 
hierarchical clustering: i) the main contributions to the growth of cluster-sized DM halos is given 
by the many lesser, closely isotropic events (see fig. 1), which are described well by the model; 
ii) such events contribute the most to statistics like the (L) — T correlation, as can be seen on 
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examining the convolutions with the complete merging histories represented by eqs. (13) and 
(14). The argument reduced to the bones goes as follows: the average (L) oc ((n2/ni)^) may be 
expressed by summing the values (n2/ni)^ after a merging episode, weighted with the probability 
n of a given mass ratio (see fig. Ic), and with the mass increment AM/M. Considering a cluster 
with M K, 10^^ Mq, the rare events (11 k, 0.15) corresponding to l^M/M k, 0.5 yield (n2/ni)^ k, 1; 
instead, events with /S.M/M ^ 0.1 have 11 ^ 0.6 and yield (n2/ni)^ ^ 6 (from eq. 9 averaged over 
the merging histories), whose product makes an overwhelming contribution to the average. 

In summary, our picture envisages the combined effect over the effective time At of all shocks 
which overlap. Barring the lumps too small to overlap and yield an appreciable mass accretion 
rate, and those too warm and too few to yield relevant compressions, our physical picture is 
focused onto a nearly continuous accretion of intermediate and colder lumps overlapping within 
a sound crossing time. From these, the ICM feels a nearly steady external pressure with only 
minor fluctuations. Such pressure exerted at the cluster boundary sets the internal density via the 
connection between boundary and centre provided by hydrostatic equilibrium, and so determines 
the steady (L), possibly varying on cosmological time scales. 

Finally, we stress the efficiency of the present analytical approach in making predictions for a 
wide range of cosmologies /cosmogonies. These are easily spanned in terms of the the convolutions 
(13) and (14) over the merging histories to yield the dependence of the amplitude, shape and 
scatter of the predicted correlations L — T, Rx — T on the cosmology as given in fig. 10 and fitted 
with eqs. 15 and 16. Such histories are dominated by those merging events between very unequal 
clumps (with Ti ^ T) which occur close to the observation time (Lacey Sz Cole 1993), as shown 
in fig. lb. Though the merging rate does depend on Qg, the sum over time of the merging events 
is weakly dependent on it, and so do the average luminosity and the dispersion (see fig. 10). 

The advances attained by the PE model are as follows. The free parameters (the central 
density Uc, (3 and 7) of the previous hydrostatic models are now computed or constrained. 
The boundary conditions that yield Uc and (3 are derived from an approximate rendition of the 
hydrodynamics, and are related to the DM dynamics. The stochastic character of the latter implies 
variance in the merging histories even at given T, and this is enhanced by the dependence of the 
emission to yield the intrinsic scatter expected in L. The resulting, narrowly constrained model 
predicts temperatures declining outwards, and - in scaling down from rich clusters to groups - 
smaller (3 and shallower brightness profiles, decreasing baryonic content on average, and the L — T 
relation bending down strongly on approaching the preheating threshold kTi ~ 0.5 keV. 

Actually, any reasonable spread in such threshold as discussed by CMT98 implies for 
groups an increased luminosity dispersion AL/L > 25% {2a) along with a considerable scatter 
A/// ~ ATi/Ti in the baryonic fraction, apart from the larger uncertainties affecting the group 
observations. We have deliberately chosen to keep the model simple and to implement here neither 
such spread nor the 2;-dependence of Ti which is expected for 2; ;^ 1, corresponding to the star 
formation rates at such early z. We plan to expand on such issues while the high-z data are 
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dr awing near. 
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FIGURE CAPTIONS 

Fig. 1. — Top panel: a Monte Carlo realization illustrates the merging history of a DM halo with 
final mass 10^^ Mq. 

Middle panel: for the same halo it is shown the probability distribution of progenitors with mass 
M' at the redshifts z = 0.03 and z = 0.5 that end up in a mass M at z = 0. The circles {z = 0.03) 
and the triangles {z = 0.5) represent the results from the Monte Carlo simulations, while the lines 

show the analytical results from eq. 4. 

Bottom panel: the fraction of merging events during the last 2 Gyrs involving a mass M and a 
partner with mass 3 10^^ < M' < M/10 (solid line); the lower limit arises from requiring a 
virial temperature Ty > 0.5 keV for the partner (see §2.2). The dashed line is the corresponding 
fraction for events with M' > M/2.5. Tilted, COBE-normalized spectrum of perturbations, as 
given in White et al. (1996) in the critical universe with Ho=50 km/s Mpc. 

Fig. 2. — The predicted dependence of (3 on the (emission-weighted) temperature T, see eq. 12 
and 8. DM potential as given by Navarro et al (1997), computed for the cosmological parameters 
and the perturbation power spectrum used in fig. 1. The shaded region indicates the 2— a scatter 
due to the merging histories. 

Fig. 3. — The predicted ratio of the baryonic fraction /2 at the cluster boundary to the external 
value fu (White et al. 1993; White & Fabian 1995) as a function of the temperature T for 7 = 1 
(solid line), 7 = 1.1 (dashed) and 7 = 1.3 (dotted). 

Fig. 4. — The ratio of the thermal ICP energy Et^^r = 3 / d^rn{r)T{r) to the ICP gravitational 

energy Egrav = Gmn J d^rn{r) M{< r)/r is shown as a function of 7 for the DM potential of fig. 
2 (solid line), and for the King form (dashed line) with core radius i?/10. The profiles n(r) and 
r(r) are provided by eq. (11). 

Fig. 5. — Temperature profiles for the model cluster of 10^^ Mq at z = in polytropic equilibrium; 
7 = 1 (solid line), 1.2 (dashed) and 1.66 (dotted), see eq. (11). DM potential as in fig. 2. The 
profile is smoothed out with a filter width of 100 kpc. The dashed area taken from Markevitch et 
al. (1997) summarizes the observations of 30 clusters. 

Fig. 6. — Upper panel: The predicted ICP density profile n(r) for a rich clusters with kT = 8 
keV (upper line), and for a poor cluster with kT = 2 kcV (lower line). Radii are normalized to the 
virial radius R. The profiles are computed using the DM potential as in fig. 2, and for the ICP 
the polytropic index 7 = 1.2. Bottom panel: to illustrate the variations produced by the use of the 
King potential as in fig. 4, with 7 = 1.1 (dotted line); and of the Moore et al. 1997 potential, with 
7 = 1.3 (dashed line). 

Fig. 7. — The surface brightness profile S(r) from the PE model is compared with the data for 
two oppositely extreme clusters. Top panel: A539, rather sparse and with low T, at 2; = 0.026 
(David et al. 1996, De Simone private communication). Bottom panel: A2390, relaxed and hot, at 
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z = 0.23 (Bohringer et al. 1998). No attempt has been made at excluding emissions from cooling 
flows. 

Fig. 8. — The predicted mass- temperature relation for Ti = 0.8, 0.5, 0.3, keV, from top to 
bottom. 

Fig. 9. — The average L-T correlation with its 2a dispersion (shaded region) is shown for the PE 
model with the tilted CDM cosmogony of fig. 1. The value Hq = 50 km/s Mpc is assumed for 
the Hubble constant. The luminosities are corrected to 0.3 solar metallicity using the standard 
Raymond Smith code. Group data from Ponman et al. (1996, solid squares); cluster data from 
Markevitch (1998, open triangles). Here Ti = 0.5 keV with no dispersion; the effect of the latter is 
shown by fig. 2 of CMT98. 

Fig. 10. — Left panel: the dependence of the average luminosity L oc {{n2 / ni)"^) on f2o. 

Right panel: the dependence of the dispersion AL on Q,o. In both panels CDM cosmogonies are 

assumed, see Liddle et al. (1996). 

Fig. 11. — The correlation of the radius Rx (see text) with T (emission weighted) is plotted for 
= 1. The 2(j dispersion is shown by the shaded region; the cosmological parameters and the 
perturbation power spectrum are as in fig. 2. The data are from Mohr &; Evrard (1997). 
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